##7b.1.

massaedados <- read.table("http://ecologia.ib.usp.br/bie5782/lib/exe/fetch.php?media=dados:babies.txt", sep = " ", header=T)

massaedados <- massaedados[massaedados$gestation != 999, ]

massaedados1 <- lm(bwt ~ gestation, massaedados)

sey <- sqrt((1/length(massaedados$bwt)+((massaedados$gestation-mean(massaedados$gestation))^2)/sum((massaedados$gestation-mean(massaedados$gestation))^2))*var(massaedados$bwt))

sey <- sey*qt(0.92,length(massaedados$bwt)-2)

posit  <- predict.lm(massaedados1, massaedados) + sey

negat <- predict.lm(massaedados1, massaedados) - sey

plot(massaedados$bwt ~ massaedados$gestation, xlab = "tempo de gestao", ylab = "peso")

abline(massaedados1, col="blue", lwd=1)

lines( massaedados$gestation, posit, type="l", col="yellow")

lines( massaedados$gestation,  negat, type="l",	col="red")


##7b.2.

init.h = c(600, 700, 800, 950, 1100, 1300, 1500)
h.d = c(253, 337, 395, 451, 495, 534, 573)

plot(h.d~init.h)
modeloa <- lm(h.d~init.h)
modelob <- update(modeloa,.~. + I(init.h^2))

anova(modeloa,modelob)

modeloc <- update(modelob,.~. + I(init.h^3))

anova(modelob, modeloc)

#O modelo mais complexo explica melhor o experimento,j que p valor foi significativo e mais baixo.


#7b.3.

massaedados <- read.table("http://ecologia.ib.usp.br/bie5782/lib/exe/fetch.php?media=dados:babies.txt", sep = " ", header=T)

massaedados.2 <- massaedados[massaedados$gestation<999,]
massaedados.3 <- massaedados.2[massaedados.2$weight<999,]
massaedados.4 <- massaedados.3[massaedados.3$bwt<999,]
massaedados.5 <- massaedados.4[massaedados.4$age<99,]
massaedados.6 <- massaedados.5[massaedados.5$height<99,]
massaedados.7 <- massaedados.6[massaedados.6$parity<9,]
massaedados.8 <- massaedados.7[massaedados.7$smoke<9,]

lm.cheio <- lm(bwt~gestation*parity*age*height*weight*smoke, data=massaedados.8)

tes0 <- lm(bwt~gestation*weight*smoke, data=massaedados.8)

summary(tes0)

tes1<-lm(bwt~gestation*weight*age, data=massaedados.8)

summary(tes1)

tes2<-lm(bwt~gestation*weight*parity, data=massaedados.8)

summary(tes2)

tes3<-lm(bwt~gestation*weight*height, data=massaedados.8)

summary(tes3)

tes00<-lm(bwt~gestation*weight, data=massaedados.8)

summary(tes00)

#o melhor modelo  com tempo de gestao e peso de me
